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ABSTRACT 


A detailed analvsis of the horizontal boundary laver structure of the warm front of 
an open ocean explosive cyclone in Intensive Observation Period (IOP) 2 of the Exper- 
iment on Rapidly Intensifving Cyclones in the Atlantic (ERICA) is conducted. Data for 
this study consists of aircraft data averaged over one minute supplimented by satellite 
and drifting buov observations. Analysis of surface winds and fluxes was done using the 
Brown-Liu Marine PBL model. Results show a PBL which differs from that found in 
typical cyclones, with large latent heat fluxes south of the warm front and with relatively 
weak sensible heat fluxes about the warm front. Boundary laver stratification was stable 
north of the warm front and unstable south of the warm front. A mechanism for moist 
frontogenesis is proposed whereby the destabilizing effects of the Jatent heat flux en- 
hances frictional convergence along the warm front. These fluxes warm and moisten the 
cvclone’s warm sector, enhancing unstable convection along the warm front and thereby 
enhancing the vertical motion. This enhanced vertical motion would strengthen the 


geostrophic deformation of the @, gradient and potentially enhance cvclogenesis. 
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I. INTRODUCTION. 


Maritime explosive cyclogenesis poses a problem for mariner and meteorologist 
alike with its unpredictable sudden severe weather. Explosive cyclogenesis is defined as 
a central pressure drop of 24 mb in 24 hours at 60°N, geostrophically adjusted for other 
=sing (Sanders and Gvakum 1980). 
sin 00 : 

Forecasting explosive cyclogenesis is difficult due to incomplete understanding of 


latitudes by multiplying the deepening rate by ( 


the triggering mechanisms that distinguish it from non-explosive development. It has 
been shown that both upper Jevel and boundary laver dynamics contribute to explosive 
development but their relative importance and the extent to which they interact is un- 
clear. It is most probable that explosive cyclogenesis has no single trigger but is rather 
the result of one of many possible combinations of forcing mechanisms. The upper level 
influences on explosive development are best known and include enhanced cyclonic 
vorticity advection due to subsidence induced vortex tube stretching (Boyle and Bosart 
1986) and enhanced frontal vertical circulation due to jet streak interaction with surface 
baroclinicity (Uccellini et.al. 1984). 

The role of boundary layer processes in explosive development is less understood, 
but preliminary investigations have focused on surface fluxes of heat and moisture as 
forcing mechanisms. Sanders and Gvakum (1980) showed that explosive cyclones were 
Winter events with a strong preference for development over the warm waters of the 
Kuroshio current east of Japan and the Gulf Stream current east of the United States. 
This relationship suggests that heat and moisture fluxes resulting from air-sea temper- 
ature differences are a factor in explosive development. Observational evidence for the 
influence of surface heat fluxes on explosive cyclogenesis was provided by Davis and 
Emanual (1988) in their investigation on the propensity of the Limited-area Fine-mesh 
Model (LFM) to underpredict rapid maritime cyclogenesis. Their hypothesis was that 
the model failed to account for warming due to surface heat and moisture fluxes. Twenty 
nine cases of explosive storms off the east coast of the United States were examined and 
in all cases the difference between predicted and actual development could be associated 
with the occurrence of warming due to surface fluxes. A concurrent climatological study 


of the same region from 1963-1977 showed a strong correlation between regions of large 


sea-surface fluxes and regions of explosive cvclogenesis. 


While emphasizing atmospheric warming due to surface fluxes, Davis and Emanuel 


(1988) did not explicitly define the flux contribution to cyclogenesis. Their contention 
was that surface heating acts to raise 1000-500 mb thickness towards a potential thick- 
ness defined by the sea-surface temperature. To achieve this thickness increase the sur- 
face fluxes must directly feed the ascending air of a cyclone. Other studies suggest 
mechanisms other than the fluxes are responsible for this atmospheric warming. Kuo 
and Reed (1988), for example, concluded that surface fluxes played a minimal role in one 
case of explosive cvclogenesis. Although they found that condensational heating asso- 
ciated with areas of heavy precipitation led to rapid development through enhanced la- 
tent heat release, the surface heat and moisture fluxes contributed verv little to this 
latent heat release. The differing impacts of the fluxes may be due to details in the dis- 
tribution of surface fluxes and boundary laver processes in developing cyclones. 

The distribution of surface fluxes in west coast non-explosive mid-latitude storms 
was examined bv Fleagle and Nuss (1985) as part of the Storm Transfer and Response 
Experiment (STREX), Surface fluxes were found to have relatively large positive values 
in the cold sector behind the cold front but small or negative values in a 200 km region 
of the warm sector in advance of the cold front. This distribution of surface fluxes would 
act to oppose cyclogenesis by reducing the baroclinicity of the cyclone through the 
warming of the cold air behind the cold front. 

This negative effect of the fluxes on cyclogenesis differs substantially from Bosart 
and Lin’s (1984) examination of the 1979 President’s Day storm, which suggested a more 
positive role of surface fluxes in explosive cyclogenesis. In this particular case, a cyclone 
developed explosively after moving offshore over the Gulf Stream in the vicinity of Cape 
Hatteras. Thev concluded that surface sensible heat fluxes alsead of the surface low acted 
to pre-condition the environment by enhancing the low-level baroclinic zone. Numerical 
model studies by Uccellini et al (1987) confirm the significant role of fluxes in fostering 
explosive development of this cyclone. The distribution of surface heating and 
moistening in this case increased low-level frontogenesis as well as directly feeding the 
regions of latent heat release in the mid-troposphere. 

Nuss (1989) examined the relationship between flux distribution and cyclogenesis 
through a series of numerical simulations using identical initial atmospheric conditions 
but different surface fluxes and sea-surface temperature distributions. The control case 
for this experiment used a zonal SST distribution with upward heat and moisture fluxes 
to the west of the surface low and downward fluxes to the east. which is similar to flux 


distributions of west coast systems. A simulation using a sinusoidal SST distribution 


to 


was then run which produced additional upward fluxes to the northeast of the surface 
iow, similar to the Presidents’ Day storm. A comparison of simulation results after 24 
hours showed that the second case deepened more rapidly than the control (1.65 mb 
h-! vice 1.47 mb h-'). The control case showed large upward fluxes in the cold sector 
behind the cold front and near zero fluxes in the warm sector in front of the cold front. 
In contrast. the second case had large upward fluxes not only in the cold sector but also 
to the northeast of the surface low in advance of the warm front. Additionally. a small 
upward flux was found throughout the warm sector. The increased deepening rate for 
the second case was attributed to the destabilizing effects of upward fluxes of heat and 
moisture on the boundary iayer which resulted in increased downward momentum fluxes 
to the north of the warm front. This in turn enhanced the associated frontal ageostrophic 
secondary circulation, leading to increased frontal convergence and vertical circulations. 
Cyclogenesis was then enhanced through an associated increase in the release of latent 
heat. 

These studies indicate that the planetary boundary laver (PBL) and frontal structure 


associated with cvclones has the following attributes: 


e A strong. stable capping inversion above a well mixed boundary laver in the cold 
region behind the cold front and in advance of the warm front: 


@® Weak or neutral stratification in the warm sector and unstable stratification in the 
cool sectors of a cyclone: 


e A narrow updraft of warm air above the surface cold and warm front position, 
Where there is no Well defined boundarv layer; 


e A frontal ageostrophic secondary circulation with the surface ageostrophic wind 
moving from the cold air to the warm. with upward motion in the warm air. 

From the lust two above items, it is clear that surface fluxes must be in the warm sector 
for them to directly feed the ascending air and warm the atmosphere. This is contrary 
to the study by Fleagle and Nuss (1985) of west coast, non-explosive storms, which 
placed the largest upward fluxes in the cold sector with small or negative surface fluxes 
in the warm sector. However, in Nuss’s (1989) simulation of the effects of surface heat 
and moisture fluxes on cyclogenesis. rapid development was associated with upward 
surface moisture fluxes in the warm sector and boundary laver stability differences across 
the warm front that apparently enhanced frontal convergence. It is not known whether 
these processes occur in actual explosive cyclones. 


Investigation into explosive cvclogenesis has been hindered by a lack of observations 


on the small scale. Most studies and their conclusions are based on synoptic scale anual- 


vsis with limited observations. This thesis will attempt to resolve the role of surface heat 
and moisture fluxes and horizontally varying boundarv layer structure by examining the 
details of the boundary layer structure of an explosive cyclone in the warm front region. 
More specifically, this thesis has the following objectives: 


I. Map the detailed horizontal boundary laver structure of an open ocean explosive 
cyclone using one minute averaged aircraft data. 7 


id 


Describe the distribution of surface heat and moisture fluxes in relation to the 
warm front. 


3. Describe the stability structure of the boundary laver across the warm front. 


4. Examine the effect that the distribution of surface heat and moisture fluxes. and 
PBL stability structure, have on the frontogenesis and ageostrophic secondary cir- 
culation of the warm front. 


Ii, DATA AND MODEL DESCRIPTION 


The Experiment on Rapidly Intensifving Cvclones over the Atlantic (ERICA) was 
devised to “determine phvsical mechanisms and processes, and their critica) spatial and 
temporal combinations, which can account for the winter time phenomenon of 
explosively developing over-ocean atmospheric storms.” (Hadlock and Kreitzberg 1988). 
Conducted from 1 December 1988 to | March 1989, the area of interest was centered 
off the east coast of the United States in a climatologicallv favored region for explosive 
cVlogenesis extending from 4U°N to 47°N between S9°W and 69°W (Hadlock and 
Kreitzberg 1988S). The ERICA field study encompassed six explosive cyclones with in- 
tensive observation periods (LOPs) of aproximately 36 hours covering the period of ex- 
plosive development. This thesis examines a portion of the second intense observation 
period (IOP-2) from 0300 UTC to 0700 UTC 14 December, which was the initial stage 
of rapid cvclogenesis. 

For JOP-2 a variety of data were obtained. The horizontal PBL structure was de- 
termined bv a NOAA WP-3D which was tasked to document the PBL s‘ructure by flying 
a sawtoothed flight profile across the warm front. Vertical structure was to be measured 
by aircraft deploved dropwindsondes, of which only one provided good data. The air- 
craft data Were supplemented by three drifting buoys from a field of fifty drifting buoys 
at 200 km spacing which measured mean sea-level pressure, air temperature, and sea- 
surface temperature (SST) on an hourly basis. Satellite SST measurements and standard 
ship observations were also available and used in this study. 

The JOP-2 aircraft data which were as used to construct detailed horizontal ana- 
Ivses, consisted of 76 meteorlogical and flight parameters recorded at one-second inter- 
vals. The first step in processing the data for this thesis was to time average the data over 
60-second intervals in order to reduce the data set to a more managable size and remove 
very small scale variablity that could not be resolved in the horizontal analysis. Care 
was taken not to average over any data gaps or around turns where wind data mav be 
questionable. Horizontal homogenity was maintained by rejecting any measurements 
taken when the aircraft was climbing or descending at greater than 5 m s"! or if the air- 
craft significantly deviated from the desired altitude of 350 m. 


A Cressman analysis scheme was then employed to produce gridded data fields from 


these processed aircraft observations. The Cressman analvsis was done over a 500 km 


by 350 km erid with a 10 km grid spacing. This grid relative to the aircraft track is shown 


in Fig 1. Five analvsis passes were made on the data using decreasing scan radu of 20.0, 
15.0, 10.0, 7.0, and 3.0 km. This thesis treats the four-hour data acquisiuon effort as a 
‘snapshot’ of the storm and no corrections for svstem translation were made. <A cor- 
rection for system translation would entail an assuruption of quasi-steady state that ‘vas 
not justified given that the cvlone deepened 12 mb over the four hours when data were 
tuken. It is believed that the errors associated with anv such correction mav be com- 
parable to analvsis errors with no correction. While this assumption about the error, 
when no translation correction is made, mav not be true in general, several factors 
minimuze the errors in this case. 


e The warm front Was quasi-staionary during the flight and the cyclone moved east- 
ward along the warm front. Consequently, structure to the east of the the low is 
not influenced by substantial translation of the primary feature and there 1s Jittle 
blending of substantially diffe.ent data, which gives relatively small errors. To the 
west of the low. the analvsis consists of pre-cyclone and post-cvclone data and the 
error is larger due to this blending: 


e The aircraft track moved to the east with the system so that the majority of the 
data are in the warm front region. In addition, the translation of the cyclone was 
relatively short compared to the data coverage to the east of the cyclone. This fa- 
vors the analysis of the warm front which ts the primary focus of this study. 

Conseyuently, the analvsis in the warm frontal region are believed to be representative 
of the structure even With some time averaging of this structure. 

The data fields were then passed to the Brown-Liu marine planetary boundary layer 
model for calculation of surface winds and wind stress, and surface heat and moisture 
fluacs. Phis PBL model is a two-laver model with an Ekman Tavlor laver over a loga- 
rithmic surface luver. The model calculates surface fluxes and winds given an upper 
boundary condition which can be the top of the PBL or within the Ekman laver. The 
flight level winds are taken as the winds within the Ekman layer from which surface 
winds are calculated. The air temperature field is required for thermal wind effects while 
the air-sea temperature difference from air and sea temperature ficlds is needed to cal- 
culate PBL straufication and surface fluxes (Brown and Liu 1982). 

The Brown-Liu model requires surface air temperature and relative humidity to cal- 
culate the imitial surface fluxes to determine the stability and iteratively solve for the 
surface fluxes. This required that the flight level data be extrapolated down to the sur- 
face. To make this extrapolation, a well mixed PBL with constant equivalent potenual 
temperature (@,) and mixing rato (q) was assumed. The temperature (7) can then be 


determined bs the relationship 


Figure 1. JOP-2 Aircraft track relative to analysis grid. 


p 286 L 
=-@{-— ahs 
1 Of [GVO ) Lp q A 


where pis the pressure at {light altitude, [ is the latent heat of vaporization, Cp is the 
specific heat at constant pressure, q is the mixing ratio, and @, is equivalent potential 
temperature measured at the aircraft. To determine relative humidity (R) at the surface, 
the saturation Vapor pressure (¢,) was calculated using the derived surface temperature 


ficld and the following relationship between vapor pressure and temperature; 


C= O.LOT8X 10” T3585 


Saturated mixing ratio (y,) was then calculated from the saturation vapor pressure and 


surface air pressure given bv the following equation: 


Finally, the relative humidity field was calculated using the aircraft mixing ratio (q) and 
the derived saturation mixing ratio (q,) using the standard formulation for relative hu- 


nudity ; 
q 
R= qn: 


To determine surface Tuxes,the model requires that the air-surface temperature and 

the air-surface humidity differences be calculated. The necessary SST field was created 

using muluple data sources in a Cressman analvsis scheme. A first guess SST field was 

derived from a NESDIS Local-Scale Analysis Region 6 gridded SST field with a 14 km 

resolution. Because the Tiros-N series polar orbiting satellites only retrieve SST’s under 

clear conditions, the age of the observations used in the analysis are variable. The par- 

ticular analvsis used was performed at 0107 UTC 14 December 1988, but the voungest 

satellite data are at least !2 hours old. and the majority of the observations were 132 

hours old. The 14 km analysis was interpolated onto a 70 x 55 10 km grid centered on ‘ 

the 40 x 35 10 km grid used for the other data fields. To eliminate any boundary effects, 

the outer 10 grid points were discarded after the analvsis was completed. é 
Three other data sources were used in the analysis: Radiometric temperatures from 

the aircraft, and sea-surface temperatures from ships and buoys. The aircraft’s down- 

ward looking radiometer observations Were then screened to eliminate erroneous data: 


I. All readings from altitudes greater than 500 m were discarded (water vapor and or 
undetected clouds in the column below the aircraft can give a cold bias); 


te 


One-munute averages were discarded if any one second had 
a. a relative humidity greater than 88°, or 
b. a liquid water content greater than 0.3 g kg-!, or 


Z 


c. a relative humidity greater than 82% and a liquid water content greater than 
O22 kev 


ta 


Data were not discarded on the basis of aircraft turns or ascent descent rate. 


Because the satellite observations were largely 132 hours old and also because the 
ship and buoy data outside the aircraft flight path were quite sparse, ship and buoy data 
dating back to 0000 UTC 9 December were examined. Ship data were included only if 
they were within approximately 3° of the satellite field or were strongly supported by 


independent aircraft or buoy duta. Older data and especially old data near strong graJi- 


ee eae eee ee 


ents were more critically examined. Most all of the buoy data were included, except in 
cases Where redundant older data were discarded in favor of collocated younger data. 
The final analvsis was performed using a Cressman scheme with a limit of 8 degrees 
on the difference from the first guess imposed as a final check on the data quality. Scan 
radi of 10.0, 5.0, 2.5, 1.0, and 0.5 km were used on five successive passes. Given the 
combination of sparse and clustered observations, these scans most closely replicated a 


hand analvsis. 


Ili. IOP-2 SYNOPTIC SUMMARY 


As discussed in Chalfant (1989), 1OP-2 is characterized by the development of se- 
veral cvclones associated with the advance of several short wave troughs in rapid suc- 
cession. The first cyclone developed well to the south and underwent a brief period of 
moderate deepening before it weakened rapidly. While the first cyclone was beginning 
to weaken, a much more significant cyclone had begun to develop to the southwest just 
off the North Carolina/Virginia coast. This second cyclone underwent verv rapid 
cvclogenesis to become the dominant feature. The 06000 UTC 13 December 1988 (from 
here on, date and time will be given as date time UTC, e.g. 13,0600 UTC) mean sca-level 
pressure and isotherm analvsis from Sanders (1989) show the low-level structure associ- 
ated with the first cvclone (Fig. 2). The cold advection evident in Fig. 2 is associated with 
the 992 mb low to the southeast of Newfoundland at 43°N, 47°W and the 1028 mb high 
off the New England coast at d8°N, 65° W. 
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Figure 2. TOP-2 13/0600 UTC surface analysis: Sea-level pressure (solid) in 
tenths of mb where the leading 9 or 10 has been dropped and surface 


temperature (dashed) in 7 (Sanders 1989). 


By 13,1800 UTC, both surface cyclones are evident as two distinct low centers em- 
bedded in a large area of cyclonic circulation (Fig. 3) just off the East Coast. The 
northern cyclone, with a central pressure of 1004 mb is located in a region of very little 
thermal gradient and very weak warm thermal advection east of the low (Fig. 3) but is 
apparently near the north wall of the Gulf Stream where a substantial surface temper- 
ature gradient exists. Differencing the isotherm analysis at 13;06 UTC and 13/18 UTC 
reveals a region of greater than 5° C warming over 12 hours to the north and east of the 
cyclone. The maximum warming occurs just east of the northern low pressure center 
at about 37.5°N, 68°W which is near the location of the Gulf Stream. Although this 
analvsis is only approximate, this rapid Warming in a region of weak thermal advection 
and near the the maximum SST gradient is certainly suggestive of warming from the 
surface fluxes associated with the Gulf Stream. Substantial upper-level vorticity 
advection is found over this northern cyclone in synoptic scale 500 mb analyses (not 


shown). Hence, there is ample support for cvclogenesis to occur in this region. 
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Figure 3. LOP-2 13/1800 UTC surface analysis: Same as figure 2 (Sanders 1989). 


At 14,0000 UTC, this northern low at 37°N, 70°W with a central pressure of 996 
mb, began explosive development. Examining the 14,0000 UTC mean sea-level pressure 
and isotherm analysis (Fig. 4), warm advection has increased at the surface just east of 
this cyclone. However, the largest surface warming is found to the northeast (40°N, 
64°W) where little warm advection is found. Extremely rapid warming of 5° C in 6 hours 
in the region to the northeast is highly suggestive of heating due to surface fluxes. A 
more substantial thermal gradient develops as the cyclone intensifies which may be due 
in part to surface warming by the heat flux. The position of the low center coincides 
with the northern boundary of the Gulf Stream as indicated by satellite SST analvsis. 
This relationship between the surface low center and the Gulf Stream was maintained 
during most of the explosive development which yielded a cyclone track along the Gulf 
Stream (Fig. 5). It is possible that such movement would enhance cyclogenesis as the 
increased surface fluxes over the Gulf Stream increases the frontogentical forcing of the 
ageostrophic secondary circulation associated with the warm front. Detailed analvsis 


using the aircraft data will explore this question in the next section. 
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Figure 4. IOP-2 14/0000 UTC surface analysis: Same as figure 2 (Sanders 1989). 


During the next 12 hours this cyclone deepened 30 mb to become a 968 mb low at 
14,1200 UTC as shown in Fig. 6. The first cyclone to the south stopped developing and 
was rapidly absorbed into this much more intense northern cyclone. The low-level 
structure at this time is typical of many mature cyclones with well defined warm and cold 


fronts and the beginning of the occlusion process. The cyclone fills during subsequent 


time periods. 


Figure 8. JOP-2 Cyclone track: SST isotherms in °C. 


To summarize, [OP-2 was a multiple low event triggered by the offshore movement 
of two upper air short wave troughs 12 hours apart. A surface low developed cast of 
Plorida in response to the first trough. The second upper air trough triggered the for- 
mation of two new lows in a surface trough that extended north of the initial svstem. 
The northernmost svstem developed explosively while moving along the northern 
boundary of the Gulf Stream in a region that had experienced rapid warming over the 
preceeding 18 hours, while the other two systems weakened and were assimilated into 
this major cyclone. Crude warming charts constructed using the temperature change 
over successive isotherm analyses suggest that surface fluxes along the Gulf Stream po- 
tenuially contributed to this extraordinary cyclone and will be explored in detail in the 


following section. 
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Figure 6. [OP-2 14/1200 UTC surface analysis: Same as figure 2 (Sanders 1989). 
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IV. BOUNDARY LAYER ANALYSIS 


From the synoptic examination in the previous section, the pattern which emerges 
is that the cold advection associated with the preceeding cyclone produced large surface 
heat and moisture fluxes that feed the developing cyclone in IOP-2. The fluxes may have 
contributed to warm frontogenesis along the Gulf Stream where the cyclone tracked. 
The aim of this section is to examine in detail the surface flux distribution and boundary 
layer structure around the warm front during the initial phase of explosive development. 
As described in section I], the aircraft data were analyzed on a grid for this analvsis, 

The detailed surface structure at around 14 0600 UTC is shown in Fig. 7. The sur- 
face low and frontal positions are taken from a separate hand analysis of the aircraft 
observations of winds and pressure. The horizontal surface air temperature distribution 
along the warm front shows a 3°C gradient from 18°to 21°C approximately coincident 
with the SST gradient. Also apparent is the cold front extending north south along 
66.5°W. The gridded thermal gradient is not exactly correspondent to the fronta) anal- 
ysis due to the data blending resulting from svstem translation to the east and smoothing 
inherent in the objective analysis method. There is, however, a high degree of confidence 
in the computed distribution about the warm front due to lack of blending as previously 
detailed in Section I]. With this in mind, analysis will concentrate on the warm frontal 
and warm sector structure of JOP-2. 

The surface heat flux distribution around the warm front (Fig. 8) shows downward 
fluxes to the north of the front and upward fluxes to the south. The magnitude of the 
downward fluxes to the north of the warm front are generally small (less than -50 W 
mr?) with more significant upward fluxes (+100 Wm-?) to the south of the front. The 
small magnitude of the fluxes north of the front indicates that the air over the cold water 
is in near equilbrium with the sea-surface temperature. The maximum heating essen- 
tially coincides with the SST gradient, which suggests the advection of cold air over the 
Gulf Stream. The increase in fluxes towards the east suggests that the air is being 
modified as it flows in from the east. This flux distribution is opposite to the distribution 
observed by Petersson et. al. (1962) and Fleagle and Nuss (1985) where cooling occurred 


in the warm air and warming occurred in the cold air. This distribution of heating resu!ts 


in Warm frontogenesis to increase the baroclinicity of the svstem. which 1s similar to the 


Figure 7. JOP-2 surface structure: Surface temperature (solid) and SST (dashed) 
in °C analyzed from aircraft data. Latitude and longitude in one degree 


increments with a center latitude and Jongitude of 38°N and 66°W. 


flux influence on the coastal frontogenesis in the Presidents’ Day storm and is suggestive 
of some arctic ice edge polar lows. 

The moisture fluxes (Vig. 9) are considerably stronger than the heat flux, with 500 
Wm”? upward fluxes south of the warm front and relatively weaker upward fluxes to the 
north. The magnitude of the fluxes in the cyclone’s warm sector are very large and ap- 
proach values recorded in cold air outbreak studies. However, these large moisture 
fluxes are occurring in the warm sector of a developing cyclone where they will feed the 
cloud and precipitation region of the cyclone. As with the heat flux, the SST and 
moisture flux gradients are essentially coincident, which emphasizes the role of the Gulf 
Stream boundary in creating the flux distribution. In addition, the moisture fiux tends 
to increase towards the east, which suggests that the easterly winds in this region are 
bringing in relatively dry air that is modified by the moisture flux. 

The fact that these strong moisture fluxes persisted throughout the initial phase of 


explosive development implies a continuous supply of cold, dry air over the warm water. 
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Figure 8. !OP-2 surface sensible heat flux: Sensible heat flux in W m2. Grid as 


in Tig. 7. 


Such a situation is not consistant with the expected warm advection from the south 
predicted bv traditional mid-latitude cvclone theory. Examination of the surface wind 
analysis (Fig. 10) reveals strong easterly flow with a wind maximum of 26 ms"! located 
south of the warm front. These strong easterly winds are consistant with the 
frontogenesis mechanism proposed by Grotjahn and Wang (1989) whereby air from the 
cold sector of a preceding cyclone is drawn into the warm sector of a following cyclone. 
The 13.0600 UTC synoptic analysis (Fig. 2) shows such a preceding cyclone to the 
northeast of the incipient explosive cyclone. A study of the boundary layer over a larger 
region that includes the tracing of air parcel trajectories would be required to confirm 
that air parcels move from the cold sector of the previous low into this region of the 
warm front. While dus is beyond the scope of this study, the wind flow near the front 
in the small scale analyses do support this interpretation. 

As discussed above, the moisture fluxes are advantageously placed to feed the warm 
frontal updraft region and thus enhance latent heat release. As suggested by Davis and 


Eimanual (1989), surface heat and moisture fluxes may act to destabilize moist svmmet- 
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Figure 9. IOP-2 surface latent heat flux: Latent heat flux (solid) in W m-?. SS] 
(dashed) in °C. Grid as in Vig 7. 


rically neutral regions and produce deep slantwise convection in the warm frontal region. 
While no analysis of the stability was done, vertical or slantwise, both satellite and radar 
images (not shown) show strong convection along the warm front with radar returns on 
the order of 50 dbz being recorded. The warming and moistening by the fluxes in this 
region may contribute to the convective instability and clearly supply moisture to the 
clouds. A more detailed analysis of the deeper structure 1s required to quantify the 
contribution of these fluxes to the convection. 

Nuss (1989) highlighted the frontogenetical role of enhanced horizontally varying 
momentum transfer to the sea surface due to the destabilizing eflects of surface fluxes 
in the boundary laver. While conclusive dropsonde data is not available for the cyclone’s 
warm sector, some idea of boundary laver stability can be inferred from comparing flight 
level winds and surface winds. South of the warm front, surface wind speeds are equal 
to flight level winds. This indicates a deep, well mixed boundary layer with flight level 
winds being mixed down to the surface, which agrees with the upward heat fluxes in this 


region. In contrast, surface winds north of the warm front are only one half the strength 
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Figure 10. TOP-2 surface winds: Isotachs in Sims! intervals. Grid as in Tig 7. 


of fight level winds, indicating a more stable boundary layer, also in agreement with the 
downward fluxes in this region. The less stable boundary laver in the cvclone s warm 
sector would enhance the transfer of momentum to the surface as suggested by Nuss 
(1989) but this increased momentum transfer is on the opposite side of the warm front 
as found in the numerical simulations. However, the effect is similar in that the resultant 
surface winds are more convergent than those at the flight level of the aircraft. This 
frictional convergence is frontogenetical which is evident from the observed structure. 

‘Phe boundary laver structure that emerges from this analysis near the warm front 
differs from that described in west coast cyclones by Fleagie et al. (1988). The unstable, 
well-mixed boundary layer occurs in the warm sector and is apparently maintained by 
the casterly flow of drv, cold air. A more stable boundary laver occurs north of the warm 
front where the air has apparently reached near equilibrium with the sea surface. This 
structure favors frontogenesis through direct diabatic heating from the sea surface and 
through its influence on the horizontally varying friction. 

To further highlight the potential impact of this boundary laver structure on the 


dvnaniucs of the cvclone, a qualitative examination of the Q-vector form of the omega 
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equation. as given by Hoskins et al. (1978) was made. As the Q-vectors will be used 


descriptively and no vectors are actually calcilated, only a brief discussion of Q-vectors 


will be provided. A complete derivation of the Q-vector form of the omega equation can 


be found in Appendix A. The Q-vector form of the omega equation is written as . 
2 4 
2 fh @ 
li er area (eae 2 ‘ 
Cp 
where o = —f,y —— is the static stability parameter and 
¢ 


= eve eV. 
Q = (0).0.) = ( Ve, ~; =*. v0 


This form of the omega equation ignores the effects of diabatic processes in forcing 
vertical mouon. If the atmosphere is saturated, then one can directly include the effects 
of latent heat release on the vertical motion by considering the forcing in terms of moist 
baroclinic dynamics. Therefore, a moist Q-vector can be defined by assuming that in a 7 
saturated atmosphere @, is equivalent to 6. Substituting @, for @ in the Q-vector omega 


equation results in : 


Q = (2.0) = Gara ae v0), 
= 


which gives the geostrophic frontogentical forcing of vertical motion in a saturated 
enviroment. A more complete descripuon of the moist frontogenesis can be obtained by 
examining the two dimensional frontogenesis equation. The moist frontogenesis 


equauion 1s 


c ( ce, | = C ( dé, | Cue CA, Cv, FG, fin €8, 
eS roa ene ec CY ueK a 


The first term represents non-conservative moist processes and is zero except at the sea 
surface Where surface heat and moisture fluxes alter 6, due to interaction with the 
boundary. The second and third terms represent the horizontal geostrophic deformation . 
of the @, gradient. which is the moist Q-vector forcing. Although the geostrophic forcing 
is not known in this case, the total wind and @, distribution suggest moist frontogentical 


forcing by the horizontal deformation (Fig. 10). In the warm sector south of the warm 
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front, large (~500 W m~) upward surface moisture fluxes and upward surface heat fluxes 
act to increase surface @, and strengthen frontogentical forcing of the warm front. This 
Warming and moistening of the warm sector also destabilizes the boundary laver in the 
evclone’s Warm sector, Which apparently results in enhanced convection along the warm 
front. This enhanced vertical motion would then strengthen the geostrophic deformation 
as indicated by the Q-vector omega equation, which potentially enhances the 
evclogenesis. While this analvsis is rather speculative. it does suggest the possible dv- 
nanucal processes by which the surface fluxes enhance cvclogenesis in this case. It re- 


mains to be shown through a more quantitative analysis whether the surface fluxes 


contribute as significantly as suggested by this qualitative analysis. 


V. CONCLUSIONS AND RECOMMENDATIONS. 


The detailed boundary laver analysis of IOP-2 has shown that the distribution of 
surface fluxes about the warm front in an explosive cyclone provides favorable 
frontogentical forcing to increase cyclogenesis. Specifically, this study has shown: 


1. The horizontal boundary laver structure of the cyclone during its explosive phase, 
when accurately mapped on a small scale using aircraft data, has a boundarv laver 
structure that differs from that found in typical cyclones. 


2. Heat fluxes about the warm front were relatively small due to air and sea surface 
temperatures being in near equilibrium. In contrast, large upward moisture fluxes 
were present south of the warm front with weaker fluxes to the north of the front. 


oo 


Boundary laver stratification was stable to the north of the warm front and unsta- 
ble to the south of the warm front. This resulted in different momentum fluxes on 
either side of the front to produce convergence in the surface winds. 


4. The surface heat and moisture flux distribution as well as the wind field are favor- 
able for moist frontogenesis along the SST gradient. Given nearly saturated condi- 
tions, this enhanced surface frontogenesis would contribute significantly to the 
evclogenesis. 

While the analvsis points to a potentially large role for surface fluxes in explosive 
cvelogenesis, the Jimitations inherent in this single case study suggest the following ree- 
ommendauons for further research: 


1. This analysis of the PBL in a limited region should be extended to a larger domain 
and include a trajectory analysis. 


2. A complete 3-D analysis should be done to examine the dynamics that are driving 
the ageostrophic circulation. 
3. The time evolution of the cyclone and PBL structure should be determined. 


4. Extend this analysis to ERICA IOPs with intermittent explosive development. A 
comparison of the surface flux distribution and magmitude during explosive and 
non-explosive phases would highlight any possible differences in the role of the 
boundary laver processes. 


APPENDIX = A. Q-VECTOR DERIVATION 


Derivation of the Q-vector form of the omega equation begins with the quasi- 


geostrophic momentum, thermodynamic and continuity equations: 


fa = : ef 
(5+ si eels Say Ue ties 


Cw 
= - — =(0 and 
Ex raat 
. 
Cu, os 
—_ +— =0 
ex cy 


For the Q-vector derivation the following assumptions are then made: 


cf 
1. f plane a 0, 
2. frictionless F, = F, = 0, 
3. adiabatic Q = 0 and 


4. constant static stability = constant. 


C 
semen cannons are then differentiated bv on and the thermodynamic 


equation by —— or —-. For the x component equation: 
c x cy 


op (e ~ — — fVog = “ and 


Equation (1) becomes 


(1) 


(2) 


—— = 3 
x (3) 
and equation (2) becomes 
G a 66 eh 6w €O ; 
(< + r,+V)r + yee + y-PS = 0. (4) 
cl . cy cy Cy cp 
Expanding the second term of equation (3), 
on ae Cu, Cu, Cy, Cu, Gu, { Cu CV, 
—.VU,= == pe 3). (5) 
ep Cp 6x Cp cy Cp Cx cy 
and by using the thermal wind relationships 
CU, a cv 5 
— = get ; f= - eee ( 
ep cy cp cy 


and the fact that the horizontal divergence of the geostrophic wind is zero, (5) reduces 
to: 


This term will be defined as Q, and is identical to, but of opposite sign to, the second 
term of (4). (3) and (4) can then be written as: 


(6) 


a 29 Au 26 
VeVi - SS =0. 7 
+0, = 2+1eo (7) 


Since the first term of (6) and (7) are equal through the thermal wind equation, sub- 
tracting (7) from (6) results in 


In quasi-geostrophic theory, geostrophic motion tends to destroy thermal wind balance 


which ageostrophic motion acts to restore. From (8) it can be seen that Q, represents 
a restoring force on geostrophic balance through differential ageostrophic horizontal 
motion acting on horizontal winds and differential vertical motion acting on the hori- 
zontal temperature gradient. 

A completely analogous derivation for the y component produces 


CUlge a6 wo _ 
So ép ” Gp éx —2Q, and (9) 
av, 
Q; = a VO. 
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To obtain the omega equation, we differentiate (8) bs, and (9) by = and add the 


two together: 


a0 3 
C Uge cA ev) CQ, 
eae ON ge os Sa nd 
cpe x CL Cx” 
a2 = 4 a 
C Vae CFA Ee 5 EQ) 
—fo = - : = -2-— 
cpes SF Sa 2 J 


Adding vields: 


Kee ~y2vo = -29-Q. 
where 
S ar, ay 
Q = (0). 0) (ar ee. ae ra 
Further simplification results in 
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whereo = —fy on. is the static stability parameter. = 
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